M'CRnror 


CH^RT 


ERGODICITY  AND  MIXING  IN  A  FIVE-MODE  SPECTRAL  MODEL  OF  TWO- 

DIMENSIONAL  TURBULENCE  DTNSRDC-S6/009 


DAVID  W.  TAYLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER 

BMtiMda,  Itorvianil  20064-6090 


CM 

CsJ 

O) 

ERGODICITY  AND  MIXING  IN  A  FIVE-MODE  SPECTRAL  MODEL 
OF  TWO-DIMENSIONAL  TURBULENCE 


< 

I 

Q 

< 


by 


M.  Calvin  Mosher 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


COMPUTATION,  MATHEMATICS,  AND  LOGISTICS  DEPARTMENT 
RESEARCH  AND  DEVELOPMENT  REPORT 


April  1986 


DTNSRDC-86/009 


I. 


MAJOR  DTNSRDC  ORGANIZATIONAL  COMPONENTS 


T 


ir 


t 


e- 


» 


s 


t- 


I 


•44 


«f,7 


NDW  DTNSRDC  3960/43b  <Rev.  2-80) 


SECURiry  CLASSlFICAnON  OF  this  page 


REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


24  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

DTNSRDC-86/009 


3  DISTRIBUTION/ AVAILABILITY  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  IS  UNLIMITED. 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME  OF  PERFORMING  ORGANIZATION 

6b  OFFICE  SYMBOL 

David  W.  Taylor  Naval  Ship 

(If  applicable) 

R&D  Center 

Code  1843 

6f  ADDRESS  (City,  State,  and  ZiPCcade) 

Bethesda,  Maryland  20084-5000 

8a.  NAME  OF  FUNDING  /  SPONSORING 

8b  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

Naval  Sea  Systems  Command 

8r  ADDRESS  (City,  State,  and  ZIP  Code) 

Research  and  Development  Office  { 

Washington,  D.C.  20362 

_ 

;b  ADDRESS  (C/fy.  State,  and  ZIP  Code) 


'0  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

61153N 


PROJECT 

NO 


WORK  UNIT 
ACCESSION  NO 

DN506133 


II  Title  (include  Security  Classification) 

ERGODICITY  AND  MIXING  IN  A  FIVE-MODE  SPECTRAL  MODEL  OF  TWO-DIMENSIONAL  TURBULENCE 


12  Personal  author(S) 

Mosher,  M.  Calvin 


1 3a  TYPE  OF  REPORT 

Final 


13b  TIME  COVERED 
FROM  TO 


14  DATE  OF  REPORT  (Year,  Month,  Day) 

1986.  April 


COSATI  COOES _ 

GROUP  I  SUB  GROUP 


18  SUBJECT  terms  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 
Dynamical  Systems 
Er god icity /Mixing 
Turbulence 


i'9  abstract  (Continue  on  reverse  if  necessary  arid  identify  by  block  number) 


^  A  dynamical  system  derived  by  truncating  to  five  frequencies  a  Fourier  representation 
of  the  Euler  equations  for  incompressible,  twoY-dimensional  fluid  flow  is  Investigated. 

In  using  Fourier  representations  to  model  turbulence,  investigators  have  assumed  that  the 
truncated  systems  are  ergodic  and  mixing,  although  the  validity  of  this  assumption  is  an 
open  question.  The  truncated  Fourier  system  has  seven  independent  constants  of  motion 
which  define  a  three-Fdimensional  torus.  Since  a  component  of  the  solution  to  the  system 
is  quasi-periodic ,  the  flow  of  the  system  can  never  be  mixing  and,  in  particular,  cannot 
be  mixing  on  the  three-dimensional  torus  defined  by  the  constants  of  motion. 


ZO  DiSTRlBurON  '  AVAILABILITY  OF  ABSTRACT 

UNCLASSIFIED  UNLIMITED  O  SAME  AS  RPT 


.'2a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

M.  Calvin  Mosher 


Zl  ABSTRACT  SECURITY  CLASSIFICATION 

□  otic  users  UNCi-ASSI  FI  ED 


ZZh  IElLPhlDNF  (Include  Area  Cod*)  22c  OFFICE  SYMBOL 

(202)  227-1932  Code  1843 


DO  FORM  1473,84mar 


83  APR  edition  may  be  usvci  until  ••.njusted 
All  other  editions  are  ohsnii'te 


OF  this  page 

UNCLASSIFIED 


The  Square  in  s*,t'  Space  Homeomorphic  to  R  with  Sides  Labeled 
by  Equations  Defining  Corresponding  Sides  of  R . 

Four  Copies  of  the  Box  R  x[0,2n]  Homeomorphic  to  a 

sq 

Three-Dimensional  Torus . 

Trajectory  ( ■—  )  of  System  Projected  onto  Region  R 

with  Boundary  ( —  ■  )  . 


ABSTRACT 


A  dynamical  system  derived  by  trxincating  to  five  fre¬ 
quencies  a  Fourier  representation  of  the  Euler  equations 
for  incompressible,  two-dimensional  fluid  flow  is  investi¬ 
gated.  In  using  Fourier  representations  to  model  tur¬ 
bulence,  investigators  have  assumed  that  the  truncated  sys¬ 
tems  are  ergodic  and  mixing,  although  the  validity  of  this 
assumption  is  an  open  question.  The  truncated  Fourier  sys¬ 
tem  has  seven  independent  constants  of  motion  which  define 
a  three-dimensional  torus.  Since  a  component  of  the  solu¬ 
tion  to  the  system  is  quasi-periodic,  the  flow  of  the  sys¬ 
tem  can  never  be  mixing  and,  in  particular,  cannot  be  mix¬ 
ing  on  the  three-dimensional  torus  defined  by  the  constants 
of  motion. 
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INTRODUCTION 

Truncated  spectral  representations  of  the  Euler  equations  have  been  used  in 

many  theoretical  and  numerical  investigations  of  inviscid,  incompressible,  two- 

1# 

dimensional  fluid  flow,  especially,  turbulent  flow.  The  representations  are 

measure-preserving  dynamical  systems  which  conserve  energy  and  total  vorticity  or 

enstrophy.  In  analogy  with  statistical  mechanics,  the  Fourier  models  are  assumed  to 
##  4 

be  ergodic  and  mixing  on  the  set  defined  by  the  conserved  energy  and  enstrophy, 

although  the  validity  of  this  assumption  is  an  open  question.  We  will  investigate  a 

5  6  7 

truncated  system  of  five  frequencies  investigated  by  Kraichnan,^  Hald,  and  Glaz' 
and,  in  particular,  the  question  of  whether  the  system  is  mixing. 

The  spectral  models  are  derived  by  representing  the  vorticity  of  an  inviscid, 
incompressible,  two-dimensional  flow  field  as  a  Fourier  series  with  complex  coeffi¬ 
cients  ujj^,  plugging  this  representation  into  the  vorticity  equation,  and  truncating 


^A  complete  listing  of  references  is  given  on  page  21. 

^Definitions  of  ergodicity,  mixing,  invariant  set,  and  quasi-periodic  can  be 
found  in  References  2  and  3. 
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the  resulting  infinite  system  of  ordinary  differential  equations  after  a  finite 

g 

number  of  frequencies.  For  each  frequency  k  the  resulting  equation  is 

du). /dt  *  0.5  1  (q  ^  -  p  Ip.qlu)  o) 

*  p+q=k  P  ^ 

2  2  2 

where  p  *  (p^,P2).  q  “  (qi,q2^  P  “  Pi'^Pa*  “  Pl^2"P2‘^1’ 

derivative  with  respect  to  time  t.  It  is  assumed  that  u)_.  is  the  complex  conjugate 

*  -K 

of  lUj^,  i.e.,  u)_j^  =  ijUj^.  The  summation  p+q  =  k  is  extended  over  only  a  finite  number 
of  frequencies.  These  systems  have  at  least  two  constants  of  motion,  namely,  the 
energy  E  and  the  squared  vorticity  or  enstrophy  Q 

-2  *  * 

E  »  0.5  Ip  uju),  n  =  Iu)U) 
p  P  P  p  P  P 

Here  the  summation  is  taken  over  all  the  frequencies  in  the  truncation.  Using  a 
truncation  which  retains  the  five  frequencies  k^  =  (1,1),  ^2  ~  “  (3,0), 

k^  =  (2,-1),  and  k^  =  (1,-1),  we  obtain  the  system 

» 

diu^/dt  *  Suj^u)^  (la) 

diu^/dt  =  (lb) 

dui^/dt  *  27(u)^u)j^  -  w^uj^)  (1c) 

dujj^/dt  =  -350)2^*  (Id) 

duj_/dt  =  -8u)_u)_  (1e) 

The  system  has  these  constants  of  motion:^ 

*  35r^  +  8r^  (2a) 


Y.,  *  35rf  +  8r^ 


'^J  fWJ  W  W  !PJ  IT.’  ly.  Al  '^J  WJ  PJ  MU  9*1  ^J.W  I  pi.f 


"y-  «  asr^  +  27r^  +  27rJ 


(2c) 


‘/^i 

ro 


Il  •  w^iu*  +  U)*U)5 

(2d) 

#  ^  » 

OvC 

Ig  «  35U)^UJ^  -  SiUgU)^ 

(2e) 

«  *  ♦. 

'-->v 

»  Im(u)^u>^u)^  +  0)20)211)5) 

(2f) 

Here  r^  *  [lo^l  for  i»l,...,5.  The  energy  E  *  enstro- 

phy  n  *  (Y^+Yg'^YjJ/SS.  However,  Y^.  Y2«  Y^*  I2’  ^3  clear  physical 

interpretation. 

Hald^  discovered  the  constants  of  motion  of  Equations  (2).  Since  the  system 
has  extra  constants  of  motion  in  addition  to  E  and  fl,  the  flow  of  the  system  cannot 
be  ergodic  or  mixing  on  the  set  defined  by  E  and  0.  However,  the  question  remains 
of  whether  the  system  is  ergodic  or  mixing  on  the  set  S  defined  by  the  new  con¬ 
stants  of  motion,  i.e.,  the  set  of  points  u)^,  w^,  lu^,  lu^^,  lu^  that  satisfy  Equations 
(2).  Also,  the  topology  of  S  is  of  interest. 

We  have  obtained  the  following  results.  The  eight  constants  of  motion 
discovered  by  Hald  are  not  independent.  A  parameterization  of  the  set  S  defined  by 
the  remaining  seven  constants  of  motion  shows  that  S  is  a  three-dimensional  torus. 
That  S  is  a  three-dimensional  manifold  implies  the  seven  constants  of  motion  are 
independent.  The  component  of  the  solution  to  the  system  can  be  solved  expli¬ 
citly  using  a  Jacobian  elliptic  function.  Since  the  u)^  component  is  quasi- 
per iodic,  the  flow  of  the  system  can  never  be  mixing  and,  in  particular,  cannot  be 
mixing  on  S.  The  results  of  a  numerical  test  are  consistent  with  the  parameteriza¬ 
tion  of  S. 

In  the  following  sections  of  this  report  we  will  first  paramterize  the  set  S 
and  then  derive  an  explicit  solution  of  the  ^^(t)  component  of  Equations  (1).  Mix¬ 
ing  will  be  discussed  next.  Finally,  the  results  of  the  numerical  test  will  be 
described. 


PARAMETERIZATION  OF  THE  INVARIANT  SET  S 
The  parameterization  of  the  set  S  defined  by  the  constants  of  motion  is  simi- 

9  i  _ _ 

lar  to  that  of  the  four-mode  case.  Letting  u).  =  r.e  for  and  i  ■  V"1f 

sJ  hJ 

we  will  show  that  the  projection  of  S  onto  r^,  r^  space  is  a  set  R  homeomorphic  to 
a  rectangle.  The  variables  r^,  r^,  r^,  8^,  8^,  8^,  and  8^  can  be  parameterized  by 
82  and  the  r^,  r^^  which  lie  in  R.  We  will  first  determine  the  set  R*  which  will  be 
the  projection  of  the  set  defined  by  Equations  (2a)  through  (2e)  onto  r2f  r^  space. 
We  will  find  that  these  equations  are  not  independent,  but  the  following  relation 
holds 


^1^2 


28011/, ‘ 


I  I 


=  0 


Therefore,  only  seven  of  the  eight  Equations  (2a)  through  (2f)  can  be  independent. 
Next,  we  will  show  that  Equation  (2f),  which  is  used  to  parameterize  8^,  defines  an 
anni.ius.  The  intersection  of  R'  and  the  annulus  will  be  R. 

i.etting  s  =  r2  and  t  =  r^^,  we  rewrite  Equations  (2d)  and  (2e)  as  follows; 

iOp-e.-ijj.+n) 


,  i(8.-82-4'J 

-1  +  r^sll/  e  =  tr^ll/i  'e 


5'‘1 


,  i(8  -8  -til.)  iO.-e.-il)  ) 

-1  >  35r^r5jl2l  ^  2  ,  2  4  2 

We  can  interpret  these  equations  by  using  triangles  in  the  complex  plane  with 
i  =  .  Let 


(3a) 

(3b) 


1 


+  iy^  be  the  value  of  the  right-hand  side  of  Equation  (3a), 


Then,  s,  t,  r^,  and  r^  satisfy  Equation  (3a)  if  and  only  if  the  triangle  ABC  in 
1  exists.  We  will  call  this  the  first  triangle.  Here  AB  is  the  complex 


r.'.rrjyer  r^sll/ 
CAB  -  arg(z.|-fl 


i(8  -8  -41/ 
e  ,  CB  is  tr^l/l 


-1  i(05-0^-4J7n) 


,  the  angle 

=  0,-82-4>,|.  and  arg(z/  =  Similarly,  let  Z2  = 


1 


be  the  value  of  the  right-hand  side  of  Equation  (3b).  Then,  s,  t,  r,,  and  r_ 

1  0 

satisfy  Equation  (3b)  if  and  only  if  the  triangle  DEF  in  Figure  2  exists.  We  will 
call  this  the  second  triangle.  Here  DE  is  the  complex  number 
,  i{0^-e  -4)2)  ..  1(02-04-4)2) 

35r^rg!l2l  e  ,  FE  is  8st|l2l  e  ,  the  angle 

FDE  *  arg(z_*l)  =  9  -8_-4»-,  and  argiz^)  =  0T-0.-4)- 


■:o 

k<a 


;.v 

I 


i 


& 


.-A 


The  equations  for  the  angles  obtained  from  the  triangle  representations  of 


Equations  (3).  i.e., 

0^  -  Sg  *  arg(2^  +  1)  +  (4a) 

0^  -  0^^  =  arg(z^)  +  -  n  (4b) 

0^  -  0^  «  arg(22+l)  +  (^c) 

02  -  0^  =  arg(22)  +  4i2 

can  be  reduced  to  the  following  four  equations: 

0^  =  02  arg(2^  +  1)  -  arg(22+1)  +  liii  ‘  >1>2  (5a) 

0j^  =  @2  -  arg(z2)  - 

0^  *  02  +  arg(2^  +  1)  +  4>^  (5c) 

arg(2^)  -  arg(2^+1)  +  arg(22+l)  -  arg(22)  '  n  (5d) 


Hence,  angles  0^,  0^,  and  0^  can  be  parameterized  by  ©2,  z^»  and  Z2.  We  will  show 
that  2^  and  22  can  in  turn  be  parameterized  by  s  and  t.  First,  however,  we  need  to 
examine  Equation  (5d). 

As  seen  in  Figures  1  and  2, 

arg(z^)  »  arg(z^-*'l)  +  angle  ABC  (6a) 


arglz^)  =  arg(z2-^l)  +  angle  DEF  (6b) 

Using  Equations  (6),  we  convert  Equation  (5d)  to  the  form 


5 


coslangle  ABC)  +  cos(angle  DEF)  *  0 


(7 


Using  the  law  of  cosines  on  the  triangles  in  Figures  1  and  2,  we  obtain 

cos(angle  ABC)  *  (r^  s^  +  r^  t^  -  |I^ {^)/(2r^r^st)  (8a 

cos(angle  DEF)  =  (35^  r^  r^  -*■  6*45^  t^  -  [ 1  ^) /(560r ^r^st)  (8b 

Using  Equations  (8)  with  r^  and  r^  expressed  in  terms  of  s  and  t  by  Equations  (2a) 
and  (2b),  we  transform  Equation  (7)  to  the  form 

Y^Y2  -  28011^1^  -  =  0  (9 

Her.ce,  only  seven  of  the  eight  Equations  (2a)  through  (2f)  can  be  independent, 
Equar.Mn  (9)  can  in  turn  be  transformed  into  the  following  three  equivalent  equa- 
tiCif.-s  uhich  will  be  useful  in  parameterizing  z^  and  z^  in  terms  of  s  and  t: 


-8V2  + 

4ll2l’^(YiV2  + 

(10a 

4ll2l'^(YiY2  * 

Ii^I^IYt  =  n20YT|iTl^|i2r^ 

(10b 

,  1  -2 

,  2  2  2  , 44 ,  ,  -2 

( 10c 

>1^2  • 

lljl  (YiY2  ^ 

To  paramietenze  z^  and  z^  in  terms  of  s  and  t,  we  use  the  following  equations 
obtained  from  the  triangle  representations  of  Equations  (2d)  and  (2e)  shown  in  Fig 


+  y-  «  35- 


did) 


Expressing  and  in  terms  of  s  and  t  by  Equations  (2a)  and  (2b)  and  using  Equa¬ 
tions  (9)  and  (10),  we  can  solve  Equations  (n)  for  x^,  x^,  y^,  and  y^  in  terms  of 
s  and  t: 


x^  *  ll^r^(y^s^  -  -  35lI^l^)/70 


(12a) 


(12b) 


y^  =  -70'^]I^1'NyiS^  +  +  2ll2|st-35ll^  1^)- 


•(Y,s^  +  Y.t^  -  2ll  'st-35ll  '^) 


(12c) 


yl  -  ibll^rSo^lI^l^^ 


(I2d) 


2  2 

We  see  that  the  equations  y2  =  0  and  y^  =  0  are  identical.  The  expression  within 
each  pair  of  parentheses  of  Equation  (12c),  when  equated  to  zero,  represents  an 
ellipse.  The  length  of  the  minor  and  major  axes  of  the  ellipses  are,  respectively, 


P_  =  2(70)’^^  1Ii|/[Yi+V2  *  (  (7,^2)^-  112011^1^)^''^ 
=  2(70)^^^  1I^|/[Yi+Y2  -  (  (71^2)^-  112011,1^)^^^ 


(13a) 


The  graph  of  the  equation  y^  =  0  is  shown  in  Figure  3.  R'  is  the  set  of  (s,t)  for 
which  the  first  and  second  triangles  exist,  or,  equivalently,  the  set  for  which  the 
right-hand  side  of  Equation  (12c)  is  non-negative.  Since  the  right-hand  side  of 
Equation  (12c)  is  negative  at  the  origin  and  at  infinity  and  s,t  are  non-negative, 
R'  must  be  the  set  of  points  lying  within  one  ellipse  but  in  the  complement  of  the 
other  in  the  first  quadrant  of  the  s,t  plane  as  seen  in  Figure  3- 

At  this  point  the  angles  6^,  6^^,  and  6^  have  been  parameterized  by  Zy  z^.  and 
@2  using  Equations  (5)  which  were  derived  from  Equations  (2d)  and  (2e).  Since  r^, 
r-,  r.,  z,,  and  z-  were  parameterized  by  the  s,t  in  R'  using  Equations  (2a),  (2b), 


(2c),  (2d),  and  (2e),  respectively,  the  set  defined  by  Equations  (2a)  through  (2e) 
is  parameterized  by  s,t,  and  6^  with  the  (s,t)  lying  in  R'  and  no  constraint  on  the 
values  of  6^.  The  only  remaining  variable  is  6^  which  will  be  parameterized  by 
Equation  (2f). 

We  will  now  determine  the  s,t  which  satisfy  Equation  (2f).  Letting 


A  »  r^t  sin(e^+0^)  -  s  r^sin(62+6^) 


B  =  -r,t  cos(6  +8,,)  +  s  rc.cos(0.+6j- ) 
1  I  *4  0  2  0 


C  =  Ij/rj 


X  =  cos(0^) 


y  =  sin(0 


3 


wf  .  ,1,.  write  Equation  (2f)  as 


C  =  Ax  +  By 


(l^a) 

(lifb) 

(l^c) 

(l^^d) 

(Tiie) 


(15) 


Sol  ■.fif’  for  we  have 

6*  =  arccosC  (AC  +  [B^(A^  +  B^  -  C^)]^^^)/(A^  +  B^)  )  (I6a) 

e'  =  arccos{  (AC  -  [B^(A^  +  B^  -  C^)]^^^)/(A^  +  B^)  )  (16b) 

Since  the  arccos  function  is  multivalued,  it  appears  that  0^  could  have  four  solu- 

t..:.;.  However,  for  fixed  values  of  A,  B,  and  C,  Equation  (15)  defines  a  line  in 

tr.c  i.ir  tables  x  and  y.  The  simultaneous  solutions  in  x  and  y  to  Equations  (lAd), 

M4ej,  and  (15)  will  be  the  intersection  points  of  the  line  and  the  circle 
2  2 

X  y  =1.  Hence,  for  fixed  values  of  A,  E,  and  C,  Equations  (i^ld),  (I4e),  and 
(15)  can  be  solved  by  at  most  two  values  of  6^  which  must  be  determined  by  Equa¬ 
tions  (l6a)  and  (16b),  respectively.  Consequently,  6^  and  6^  are  single-valued. 
When  the  line  is  tangent  to  the  circle.  6^  has  one  solution  and  6^  equals  0,, 


i.e.,  0,  has  two  solutions  if  A  +  B  -  C  is  positive  and  one  solution  if 

2  2^  2 
A  +  B  -  C  is  zero. 

Finding  the  s,t  which  satisfy  Equation  (2f)  is  now  reduced  to  finding  the  s,t 

2  2  2  2  2 
for  which  A  +  B  -  C  is  non-negative.  Defining  u  ■  s  +  t  and 

P(u)  »  A^  +  B^  -  C^,  we  have 

P(u)  «  -8u^/35  +  (Yi+Y2>u/35  -  |I^j^  -  35|l3[^/(Y3-27u) 


We  will  show  that  the  values  of  u  satisfying  Equations  (2)  must  lie  between 

two  roots  of  P(u)  in  the  interval  P^]n(0,Y3/27) .  Let  Q{u)  »  {Y3-27u)P(u) . 

Since  0(u)  is  a  cubic  polynomial  in  u,  it  has  three  roots.  Consequently,  P(u)  also 

2 

has  three  roots  because  Q(-i^/27)  *  "35 113 1  with  I3  being,  in  general,  nonzero. 

Since  P(u)  is  positive  as  u  approaches  Y2/27  from  the  right  but  negative  as  u 

approaches  +®,  P{u)  has  one  real  root  greater  than  'i^/27. 

However,  the  two  remaining  roots  of  P(u)  must  lie  in  the  interval 

P^] n(0,Y3/27) .  In  order  that  Equation  (15)  be  solvable  for  nontrivial  values 

of  63,  P(u)  must  be  positive  for  some  value  of  u'  corresponding  to  a  point  (s',  t') 

in  R'.  In  particular,  u'  must  be  less  than  •i^/27  to  satisfy  Equation  (2c).  Hence, 

^  2  2 

P(u)  must  be  positive  for  some  values  of  u  in  the  interval  [P_,  P^]n(0,Y,/27). 

2  2  ^ 

Since  P(u)  is  negative  at  both  p_  and  p^  and  is  also  negative  as  u  approaches  Y3/27 

from  the  left,  P(u)  must  have  its  two  remaining  roots  u^,  u^  in 

[P?,  P^]n(0,Y3/27) .  Hence,  P(u)  is  non-negative  for  (s,t)  in  R'  which  lie  in  the 

annulus  with  inner  radius  and  outer  radius  ^u^  assuming  u^  is  less  than  u^. 

We  conclude  that  the  s,t  for  which  6-  is  defined,  i.e.,  the  s,t  for  which 

p  p  2  ^  —  “ 

A  +  B  -  C  is  non-negative,  lie  in  the  annulus  with  radii  ^u^  and  ^u^.  The  set  R 
shown  in  Figure  is  the  intersection  of  this  annulus  and  R' .  R  is  the  set  of 
points  in  s,t  space  which  we  use  to  parameterize  the  set  S. 

The  coordinates  6^,  63,  0^^,  and  6^  of  the  points  in  S  have  now  been  parameter¬ 
ized  by  the  coordinates  s,t,  and  0^  of  points  in  the  set  R  x  [0,2n].  Equations  (5) 
are  used  to  parameterize  8^,  6^^,  and  0^  and  Equations  (16)  are  used  for  63.  Equa¬ 
tions  (2a),  (2b),  and  (2c)  are  used  to  parameterize  r^,  r^,  and  r3,  respectively. 
Collecting  these  equations,  we  have: 


«  t{V,-8r)/35] 


(17a 


^2 

"3 

^5 

«2 


5 


6- 

< 


t 


4' 


S 


[(Y2-27(s^+t^))/35]^^^ 

t 

[(Y2-8t^)/35]^^^ 

©2  +  arg(2^  +  1)  +  4i^ 

®2 

©2  -  &rg(z^)  -  4I2 
©2  +  arg(2^+l)  -  arg(22+l) 
arccos^  (4>^) 

arccos2(‘('*) 

1arccos^(<l>  ) 

arccos2('1’  ) 

(AC  +  [B^(A^  + 

(AC  -  [B^(A^  +  B^  - 


-  4.1  -  4«2 

if  \(arccos^  ) 

if  X(arccos2(<t>'^) ) 
if  X(arccos^ (t  ) ) 

if  X(arccos2('1>  ) ) 
)/(A^  +  B^) 

)/(A^  +  B^) 


*  0 

=  0 

=  0 

=  0 


(17b) 

(17c) 

(17d) 

(17e) 

(17f) 

(17g) 

(17h) 

(17i) 

(17j 

(17k 

(171 

(17m 


X(0)  =  C  -  A  cos(0)  -  B  sin(0) 

Here  \(0)  =  0  is  Equation  (2f)  and  arccos^  for  i  =  1,2  are  the  two  branches  of  the 
arccos  function  taking  values  between  0  and  n  and  between  n  and  271,  respectively. 


The  variables  A,  B,  and  C  are  defined  by  Equations  {14a),  (I4b),  and  (I4c),  respec 
tively,  and  and  are  defined  by  Equations  (12). 

Equations  (17),  however,  cannot  be  used  in  their  present  form  as  a  one-to-one 
parameterization  of  S  because  0^  has  multiple  values.  In  addition,  the  signs  of  y^ 
and  y^  are  undetermined  implying  both  branches  of  the  arg  function  must  be  used  in 
Equations  (17f),  (17h),  and  (17i).  As  indicated  in  the  derivation  of  Equations 
(16),  Sj  has  two  values  for  any  <  s,  t,  Og  >  in  R  x  [0,  27i],  i.e.,  and  0j. 

Also,  for  any  <  s,  t  >  in  R  Equations  (17)  do  not  specify  the  signs  of  y^(s,t)  and 
y2(s,t)  implying  arg(z^)  and  arg(z^+1)  for  i  *  1,  2  will  each  have  two  values  for 
any  <  s,  t  >  in  R.  However,  only  the  sign  of  y^  needs  to  be  determined  because  we 
can  show  that  y^  and  yg  have  opposite  signs  for  any  s,t  in  R.  The  argument 
proceeds  as  follows.  Since  arg(z^)  =  arg(z^+1)  +  angle  ABC  as  seen  in  Figure  1,  w 
have 

arg(z^)  -  arg(z^+1)  >  0  for  y^  >  0 
and 

arg(z^)  -  arg(z^+1)  <  0  for  y^  <  0 

Exactly  the  same  inequalities  hold  with  z^,  y^,  and  angle  ABC  replaced  with  z^, 
y^,  and  angle  DEF,  respectively,  in  Figure  2.  If  y^  >  0,  the  angle  ABC  lies 
between  0  and  n  implying  by  Equation  (5d)  that 

-  arg{z2+1)  •  angle  DEF  <  0 

and,  therefore,  y2  <  0.  Similarly,  if  y^  <  0,  then  y^  >  0. 

To  obtain  a  one-to-one  parameterization  of  S,  we  can  account  for  the  two 
values  of  0^  and  the  sign  of  y^  by  reformulating  Equations  (17)  as  four  maps 
ijj^,  ...,  corresponding  to  four  cases: 

0^  =  Bjt  y-|  i  0  for  case  1  (I8a 

0-  =  0!.  V.  <  0 


for  case  2 


{18b 


03  ■  03'  'i  -  ° 


for  case  3 


(18c) 


63  *  63,  <  0  for  case  il  (18cl) 

Each  4)^  is  a  one-to-one  map  from  R  ^  [0,2n]  to  S 

4). :  R  X  [0,  2n]  ---►  S 

for  i  =  The  map  4<i  is  defined  by  Equations  (17)  with  63  ■  63  and  with 

the  arg  function  in  Equations  (17f),  (I7h),  and  (17i)  taking  values  between  0  and 

n,  i.e.,  y^  >  0.  The  map  will  also  have  values  63  defined  by  63  but  the  arg 

function  will  take  values  between  n  and  271,  i.e.,  y^  <  0.  However,  4*^  will  have 

values  63  defined  by  63  and  the  arg  function  defined  for  >  0.  The  map  4*2^  will 

have  values  63  defined  by  63  and  the  arg  function  defined  for  y^  <  0. 

To  construct  the  parameterization  4)  of  S  from  41^ ,  . . . ,  4j^  we  need  to  paste 

together  four  copies  of  R  x  [o,27t].  To  facilitate  the  construction  of  4>»  we  can 

view  R  X  [0,2/1]  as  a  box  because  R  is  homeomorphic  to  a  square.  That  is,  R  can  be 

mapped  by  a  homeomorphism  onto  a  square  R  in  a  space  with  new  coordinates  s',t' 

sq 

as  shown  in  Figure  5;  the  two  sides  of  R  determined  by  the  ellipses  in  Figure  4  are 

mapped  to  opposite  sides  of  the  square  and  the  two  sides  determined  by  the  annulus 

are  mapped  to  the  other  two  sides  of  the  square.  The  four  copies  of  R  x[o,27i] 

sq 

viewed  as  boxes  can  be  pasted  together  with  corresponding  sides  matched  as  shown  in 
Figure  6  to  form  a  composite  box  with  periodic  sides,  i.e.,  a  three-dimensional 

torus.  Each  map  4<^  will  be  defined  on  the  i^  box  in  Figure  6  for  i=l,...,4.  On 

the  surfaces  of  the  boxes  the  4*.  fof  i=l,...,4  have  values  in  common  because 

+  -  ^  2  2 
6-  =  6-  for  any  <  s,  t.  Bp  >  with  <  s,  t  >  on  the  circles  defined  by  s  +  t  =  u. 

and  s  +  t  *  U2  in  Figure  4.  Also,  the  two  branches  of  the  arg  function  (for 

y,  >  0  and  y,  <  0)  agree  for  any  <  s,  t,  6,  >  with  <  s,  t  >  on  the  ellipses  defined 

by  =  0  in  Figure  4.  Since  each  4*^  is  one-to-one  and  continuous  on  the  i  box 

and  the  4)^  have  values  in  common  only  on  the  surfaces  of  the  boxes  for  i=l,...,4, 

we  conclude  that  the  4*1 »  •••»  define  a  nomeomorphism  4^  from  the  three- 

dimensional  torus  consisting  of  four  copies  of  Rg|^x[o,27i]  pasted  together  to  the 

set  S  defined  by  the  constants  of  motion. 


EXPLICIT  SOLUTION  OF  lu^Ct) 

An  explicit  solution  of  the  lu^Ct)  component  of  Equations  (1)  will  be  derived. 
We  start  by  differentiating  Equation  (1c)  with  respect  to  time  t.  Using  Equations 
(2a),  (2b),  and  (2c)  we  obtain 

d^Uj/dt^  »  (a+Prg)^^  (19 

where  a  *  iby^  "  27(V^+Y2^  ^  *  “560.  Converted  to  polar  coordinates.  Equation 

(19)  can  be  represented  as  two  equations 


d^r^/dt^  -  pr^  -  ar^  “  r2(de2/dt)^  =  0  (20a) 

2(dr2/dt)  (dO^/dt)  +  r^  (d^P^/dt^)  =  0  (20t 

Integrating  the  second  equation  with  respect  to  t,  we  obtain 

dPj/dt  =  r^^  e  ’ 

where  c,  is  a  constant.  If  we  substitute  d6,/dt  into  the  first  equation,  multiply 

5  2 

by  dr^/dt,  integrate  and  multiply  the  result  by  r^,  the  first  equation  becomes 


2C 

(dy/dt)^  »  2Py^  +  -  Sc^y  -  ^e  ^ 


where  y  *  r^  and  c^  is  a  constant  from  the  integration.  Equation  (22)  can  be  writ 
ten  in  the  form 

(dy/dt)^  «  h^(y-X^)(y-\2)(y-\2)  (23 

2 

where  h  *  2p  and  X,,  X,,  and  X-  are  the  roots  of  the  right-hand  side.  This  is  an 

'  ^  10 
elliptic  equation.  It  has  an  explicit  solution 

r^  =  y  *  X^  +  (X2-X2)sn^(hM(t+tp) ,k)  (24 

P  2 

where  M  »  0.25(X^-X2),  k  *  (X2-X2)/(X^-X2) ,  and  sn  is  a  Jacobian  elliptic  func¬ 
tion.  Hence,  r,  is  periodic  with  period  p  «  2K/hM,  where 


Now  6^(t)  can  be  obtained  by  Integrating  Equation  (21)  with  respect  to  t  using 
Equation  (2^). 

MIXING 

We  will  show  that  the  real  and  imaginary  components  of  ^^(t)  are  quasi- 
periodic.  Letting  and  y^  be  the  real  and  imaginary  parts  of  uj^,  respectively, 
we  can  convert  Equation  (19)  into  two  second  order  equations 

d^x^/dt^  =  (a+pr^lx^  (25a) 

d^y^/dt^  =  (o+Pr^ly^  (25b) 

Letting  v^  =  x.  and  v.  =  dx./dt,  we  now  convert  (25a)  into  a  first-order  system  of 
two  equations  ^  ^ 

dv/dt  =  H(t)v  (26) 

-  2 

with  V  ■  (v  ,  V-)  and  H(t)  a  2><2  matrix.  Because  r,(t)  is  periodic  with  period  p, 

H(t)  is  periodic  with  the  same  period.  By  Floquet  theory  there  exists  a  periodic 

nonsingular  matrix  P(t)  of  period  p  and  a  constant  matrix  R  such  that  a  fundamental 

matrix  4i(t)  of  the  system  (26)  can  be  represented  as 

♦(t)  =  P(t)e'^'^ 

However,  the  eigenvalues  of  the  matrix  e  must  be  pure  imaginary  because 

r2(t)  =  making  X2(t)  =  Equation  (25a)  and  both  x^  and  dx^/dt  are 

bounded  as  seen  by  inspection  of  Equations  (Ic),  (2a),  (2b),  and  (2c).  Transform- 

R  12 

ing  e  into  Jordan  canonical  form  and  representing  the  transformed  matrix  in 

exponential  form  e  ,  we  see  that  the  new  fundamental  matrix  it)  =  P'(t)e  (P' 

being  the  transformed  matrix  P)  must  be  quasi-periodic.  Hence,  ^^(t)  is  quasi- 

periodic  as  is  y^ft)  by  the  same  argument  applied  to  Equation  (25b). 

Since  ^^(t)  is  quasi-periodic,  the  solution  to  Equations  (1)  can  never,  in 

general,  be  mixing.  In  particular,  since  is  periodic,  the  solution  can  never 

be  mixing  on  the  set  S  defined  by  the  constants  of  motion. 


NUMERICAL  TEST 


The  Gear-Hlndmarsh  ^  ordinary  differential  equation  solver  was  used  to  solve 

the  system  of  Equations  (1)  for  the  initial  point  ui,  ■  1 ,  u).  »  1 ,  u>-  *  1 ,  Uj.  »  1, 

I —  -7  '  ^  ^  ** 

and  u)^  >  v-1>  The  parameter  EPS  was  set  at  10  .  At  every  10  time  steps  the  solu¬ 

tion  was  plugged  into  Equations  (2)  to  ensure  that  the  resulting  constants  of 
motion  agreed  to  four  digits  with  those  of  the  initial  point,  namely,  Yi  ”  ^3i 
>2  •  Y3  »  89,  «  1  +  V^,  Ig  »  -8  -  35V^,  and  •  -1.  An  Apollo  DN420  com¬ 

puter  was  used  with  fifteen  digits  of  accuracy  (double  precision). 

When  the  system  was  integrated  to  t  *  1 ,  the  trajectory  of  the  system  wound 
around  the  torus  passing  repeatedly  through  the  four  boxes  in  Figure  6.  The  tra¬ 
jectory  projected  onto  s,t  space  intersects  repeatedly  the  boundary  of  R  but  never 
crosses  it  as  seen  in  Figure  7.  That  the  trajectory  appears  in  Figure  7  to  be 
tangent  to  the  boundary  of  R  while  never  crossing  it  is  evidence  that  the  deriva¬ 
tion  of  the  region  R  is  correct. 


7  7  2  2 

5^  ♦  r  *  Uj  ♦  r  =  U2 

Figure  6  -  Four  Copies  of  the  Box  Rsq>‘[0,2x]  Homeomorphic 
to  a  Three-Dimensional  Torus 
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